Sum rule for the pair correlation function 
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A sum rule has been derived for the static pair correlation function. This rule is the extension of 
the well-known equation that relates density fluctuation to compressibility. The obtained sum rule 
is applied to the Bose and Fermi ideal gases as well as BCS and Bogoliubov's models. 
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The well-known sum rule that relates density fluctua- 
tion to compressibility of the equilibrium thermodynamic 
system plays important role, since it is applicable to de- 
scription of classical as well as quantum liquids at all den- 
sities. The relation, obtained by Ornstein and Zernike as 
early as in twenties for a homogeneous system (see, e.g., 
Ref. [0 and references therein), reads 



d^rn[g{r) — 1] = 



T 



dn 



-1, 



(1) 



where g{r) is the pair distribution function, which de- 
scribes the density correlations, and T, fi, and n = N/V 
denote temperature, chemical potential, and density of 
particles, respectively. The first term in r.h.s. of Eq. (|l|) 
is proportional to the fluctuations of the number of par- 
ticles in the Gibbs grand ensemble and can be connected 
also with the thermal compressibility: T{dn/ d^)T /n — 
T{dn/ dp)T = {6N'^)/N, here the flrst equation is the 
thermodynamic one, and the second one results from 
Eq. (^ for A = N (see below). Below we shall prove 
that the following generalization is valid for a thermody- 
namic system in the equilibrium state: 



d2 [(v;t(i)V't(2)V'(2)vi(i')>"^ 

- (V3t(2)^(2))-(^t(l)^(l'))a 



where {■ ■ ■)^^ means the average over the Gibbs canon- 
ical grand ensemble in the thermodynamic limit 
1 = (ri,cri) stands for the coordinate and spin indices 
of a particle, / d2 • • • = J2a2 1 '^'^''2 • • •, and -ip denote 
the Bose or Fermi field operators. This relation is cor- 
rect even for a non-homogeneous system. In the homo- 
geneous case the local density of particles {i/;'' {2)tp{2))^ 
becomes constant, and r.h.s. of Eq. (||) depends only on 
ri — r'l , which allows us to employ the Fourier transfor- 
mation when investigating a specific model. We stress 



that the sum rule contains asymptotic limiting values of 
the correlation functions, and, hence, one can replace 
the condensate operators by c-numbers in the case of 
Bose system below the temperature of the Bose-Einstein 
condensation, as was proved by Bogoliubov The re- 
lation (^ is the generalization of Eq. (^, which is ob- 
tained from Eq. (^ by setting 1 = 1' and using the for- 
mula {ip'' {l)tp {!)) — n; here we put by definition 

5(|r2 - rij) = Eai,..(V'Ul)^U2)^(2)^(l))^V"' for the 
pair distribution function. 

The proof is rather simple. Let us consider a subsystem 
of volume V embedded in a system of volume Vq with 
which the subsystem can exchange energy and particles. 
It follows from the definition of the number operator TV = 
Jy d2ip'^ {2)ip{2) of the subsystem and the commutation 
relation [iV,^^(l')] = -ip{l') 



d2 [(V't(i)^t(2)^(2)V-(l')) 

V 

-(V't(2)V'(2))(^t(i)^(l'))] 

{^\im'){N - N)) - iij^iim')), (3) 



here the coordinates ri and r'^ in the indices 1 and 1' 
range over the subsystem only. Let us perform the ther- 
modynamic limit Vq oo, n — const; then the subsys- 
tem can be considered as the grand ensemble. Taking 
after that the limit V oo and using the general rela- 
tion for the grand ensemble {[A, TV] = is assumed) 



(AN) ~ {A){N) ^Td{A)/dii 



(4) 
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for A = V-Ul)^!!'); we arrive at Eq. (|). Note that 
performing the limit Vq ^ oo is the essential point of the 
proof, for the two-body correlation function in Eq. (^) 
can contain, apart from its asymptotic value, terms of the 
order of 1/Vb, which do not tend to zero at large distances 
between ri, v'^ and T2 and could give a finite contribution 
to the integral in Eq. (||). The identity (jj) can be verified 
directly from the expression of averages {A) = Tr {Ap) in 
the canonical grand ensemble with the density matrix p = 
iexp[-(i7 - pN)/T], here Z = Tr exp[-(7? - pN)/T] 
is the grand partition function [Q . 

It is not difficult to convince ourselves that the sum 
rule (H) is satisfied for the ideal Fermi gas or Bose one 



1 



above the critical temperature. Indeed, using the Wick's 
theorem and performing the Fourier transformation yield 
^n^^ = T9nko-/9/i — nko-, which is obviously fulfilled 
for the occupation numbers. Below the critical temper- 
ature divergences appear in Eqs. (|l|), (||) and (^). This 
is due to the fact that the chemical potential of the ideal 
Bose gas becomes zero and independent of the density 
and temperature in the thermodynamic limit; as a con- 
sequence, the derivatives with respect to the chemical 
potential go to infinity. In particular, that means that 
the fluctuations of the total number of particles are non- 
thermodynamic in the grand ensemble: {SN'^)/N oo 



in the limit V 



instead of the normal behaviour 



N, which Eq. (W) implies for A N; this is 



well-known anomaly of the ideal Bose gas Q . 

The correlation functions in the exactly solvable model 
of Bardeen, Cooper and Schrieffer can be evaluated also 
with the Wick's theorem [^, as the Hamiltonian can be 
approximated (exactly in the thermodynamic limit) by 
the quadratic form of the Fermi operators. It is not dif- 
ficult to see that in this case the sum rules (|l]) and (^) 
are violated. The violation of Eq. (|l|) was shown by J. 
Bell, who revealed parallels between this problem and 
the well-known problem of the gauge-invariant Meissner 
effect in BCS model §. 

Let us now examine the Bogoliubov model for spin- 
less bosons. At first, we select the c-number parts in 



the field operators ip' 



"d^ and tp 



z9. 



Then the correlation functions are calculated easily with 
the Wick's theorem for the operators and i9, since 
the model Hamiltonian is the quadratic form of the 
Bose operators Q. Thus, in addition to obvious rela- 
tions (i?) = (??^) — 0, we have for the triple averages 
= {-d^d^d) — 0. Performing the Fourier transfor- 
mation of Eq. (0) yields for fc = 



2nQ lim {^pi^ 



-no 



T- 



1 



(dfi/dno) 



and for fc 7^ 



(5) 



(6) 



Here '0k = (akO-k) and rik = (akflk) stand for the nor- 
mal and anomalous averages. The relation (^, which cor- 
responds to the condensate density, is satisfied exactly; 
this can be verified with the expressions of paper |^ with 
fj, = no<&(0) (here $(0) is the zero Fourier component 
of the interaction potential). At the same time, one can 
readily see that Eq. (||) is not fulfilled at zero temperature 
due to the terms IV'kP and rt^, resulting from the four- 
boson average {d^ di)) . At non-zero temperatures we 
arrive at obvious divergence in the l.h.s. of Eqs. (0) and 
(||) due to those terms, since I'f/'kP l/k^ aX k ^ 

according to the Bogoliubov's theorem [||. Note 

that the three- and four-boson averages are supposed to 
be of the next order of smallness in comparison with j-i/jkl 



and rik within the Bogoliubov theory; in particular, they 
are neglected in the Bogoliubov's Hamiltonian ||^ . From 
Eq. (|6[) one can conclude that evaluation of the four- 
boson averages is not consistent in the framework of the 
Bogoliubov model and, consequently, exceeds its accu- 
racy. It is worth noting that in calculating the pair dis- 
tribution function g{r) Bogoliubov restricted himself by 
the two-boson terms ipk and rik only 

In conclusion, the generalization (^ of the well-known 
relation (|^) is proposed. Violation of the sum rule (H) 
may indicate an anomaly in a model or point out limits 
of its validity. In particular, we show that calculation of 
the four-bosons averages is beyond the accuracy of the 
Bogoliubov model. 
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